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We present results of a microscopic density functional theory study of wedge filling transitions, at 
a right-angle wedge, in the presence of dispersion-like wall-fluid forces. Far from the corner the walls 
of the wedge show a first-order wetting transition at a temperature T^ which is progressively closer 
to the bulk critical temperature Tc as the strength of the wall forces is reduced. In addition, the 
meniscus formed near the corner undergoes a filling transition at a temperature T/ < T^, the value 
of which is found to be in excellent agreement with macroscopic predictions. We show that the filling 
transition is first-order ii it occurs far from the critical point but is continuous if T/ is close to Tc even 
though the walls still show first-order wetting behaviour. For this continuous transition the distance 
of the meniscus from the apex grows as £w ~ (Tf — T)~^^ with critical exponent fi^ ~ 0.46 ± 0.05 
in good agreement with the phenomenological effective Hamiltonian prediction. Our results suggest 
that critical filling transitions, with accompanying large scale universal interfacial fluctuation effects, 
are more generic than thought previously, and are experimentally accessible. 
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There is now direct experimental evidence for the ther- 
mal excitation of the gravity stabilised capillary-wave- 
like fluctuations at the interface between coexisting fluid 
phases [1 . Over the last few decades theory has predicted 
that such fluctuation effects are particularly important at 
certain types of interfacial phase transition such as crit- 
ical wetting [IHS]- Wetting refers to the unbinding of 
a fluid interface from a solid substrate (or another fluid 
interface) on approaching a temperature T^,, at which 
the contact angle 9 vanishes. The order of these transi- 
tions is determined by the subtle interplay between wall- 
fluid and fluid-fluid intermolecular forces and also inter- 
facial fluctuations. The original macroscopic argument 
for wetting transitions had predicted that the transition 
would be first-order, and should necessarily occur on ap- 
proaching the bulk critical temperature Tc [7j. Model 
calculations soon revealed that the location and order of 
the transition are more general than this. In particu- 
lar, Nakanishi and Fisher showed that, for systems with 
short-ranged forces, the transition should change from 
first-order to continuous, if the surface forces are weak- 
ened and Tu, approaches Tc [S]. While this has been 
fully tested in Ising model studies [U], this scenario is 
altered by the presence of long-ranged, dispersion-like 
intermolecular forces. In order to see continuous (now 
referred to as "critical" ) wetting transitions one requires 
a fine tuning of the range and strengths of the solid-fiuid 
and fluid-fluid forces [21 [10] . Consequently while there 
are many examples of first-order wetting, there are no 
unambiguous experimental examples of critical wetting 
for solid-fluid interfaces, although the transition has been 
seen in a few binary mixtures [5|. 

One way around this, which would allow one to see the 
strong influence of interfacial fluctuations on a continu- 
ous phase transition, is to consider fluid adsorption in a 
linear wedge for which there is an analogous transition 
referred to as filling [TTMlSj . This transition is far more 
common in nature than the wetting transition and was 



first studied experimentally almost 40 years ago [16] al- 
though the order of the transition was not considered. 
Recent phenomenological effective Hamiltonian models 
have predicted that fluctuation effects are enhanced com- 
pared to wetting and also that the requirements that the 
transition can be continuous are more relaxed |17H19j . 
While this has been studied extensively in the Ising model 
20 - 24J , the more realistic case of long-ranged forces has 
not been studied in detail. In this paper we present the 
results of a study of filling in the presence of dispersion 
forces, based on a microscopic classical density functional 
theory (DFT). The latter has been instrumental in devel- 
oping our understanding of inhomogeneous fluids but is 
most usually applied to systems in which the equilibrium 
density depends on only one co-ordinate [25H25] . Here we 
use a 2D DFT to study filling transitions and compare 
with the predictions of thermodynamic arguments and 
effective Hamiltonian theory. We find that close to Tc 
the filling transition is continuous even though the walls 
of the wedge themselves still exhibit first-order wetting. 
This allows us to check interfacial Hamiltonian predic- 
tions for the critical behaviour and is strongly encourag- 
ing that continuous filling transitions may be found in 
the laboratory similar to experiments on complete wedge 
fifling [ig. 

Consider a wedge geometry formed by two identical 
infinite planar walls that meet at an opening angle 2il) in 
contact with a bulk vapour at chemical potential /^, tuned 
to saturation /i — /i^^ at a temperature T < Tc- The 
wedge may be thought as being a missing link between 
a planar wall [i) = 7r/2) and a capillary-slit {if} = 0) and 
shows a phase transition, filling, which is distinct from 
wetting and capillary condensation. Far from the apex 
the thickness of the liquid wetting layer, £^, is the same 
as for a planar wall. However near the apex, the thickness 
of the meniscus can be much greater. Macroscopic argu- 
ments dictate that the wedge is completely filled above a 
filling transition temperature Tf which occurs when the 



contact angle of a liquid drop satisfies [TTI - fT^ 



(1) 



The wedge filling transition corresponds to the change 
from microscopic to macroscopic adsorption, asT ^ Tf, 
and may be first-order or continuous (critical filling) cor- 
responding to the discontinuous or continuous divergence 
of the adsorption. Because (ll]) is an exact requirement, 
the filling transition is ubiquitous in nature for all fluids 
that have form drops with a finite contact angle. 

Within classical DFT the equilibrium density pro- 
file is found by minimizing the grand potential func- 
tional n[p] = F[p] + Jdrp{r)[V{r) - p], where V^(r) is 
the external potential [3D|. We consider a right angle 
wedge {tp = 7r/4) so that the potential V{r) = V{x, z) 
is a function of Cartesians x, z > and is transla- 
tionally invariant along the wedge. Here F[p] is the 
intrinsic free energy functional of the fluid one-body 
density, p(r), which can be split into ideal and excess 
parts. Modern DFT often divides the latter into a hard- 
sphere part Fhs [p] and an attractive contribution Fa [p] = 
^ / /dridr2/o(ri)p(r2)ua(ri2) where Ua{r) is the attrac- 
tive part of the fluid-fluid potential. We take this to be a 
Lennard- Jones (LJ) potential Ua(r) = — 4e((T/r)^iJ(r— ct) 
which is truncated at r^, = 2.5 u, where a is the hard- 
sphere diameter and H{x) is the Heaviside function. For 
Fhs[p] we use Rosenfeld's fundamental theory which ac- 
curately models packing effects if the density is high close 
to the walls [31] [31] . The external potential arises from 
a uniform distribution of wall atoms, with density p^, 
which for r > a interact with the fluid atoms via the LJ 
potential 0„(r) = —4:8.^ (-) , leading to 



V{x, z) — an 
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2z4 + x^z^ + 2x'^ 
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2x^z^^/x'^ + z2 x^ 



(2) 



where a^ = —-oT^£wPw<^^- There is a hard- wall repulsion 
if X, z < a. Infinitely far from the wedge apex, the po- 
tential close to either surface recovers that of a planar 
wall e.g. y(cxD, z) — loiruij z^ . The functional ri[p] is min- 
imized numerically on an L x i grid where the lateral 
dimension of our box is L = 50cr and the grid has dis- 
cretization size 0.05 a. To mimic the bulk boundary con- 
ditions we impose p{L,z) = Ptt{z) and p{x,L) = pTr{x) 
where PTr{z) is the equilibrium profile for a planar wall- 
fluid interface with p^riL) fixed to the bulk gas density 
Pg. In our model DFT ksTc/e — 1.414 and temperature 
is expressed either in fractions of T^ or in dimensionless 
units T* = kBT/e. 

We have considered a variety of wall strengths and 
present results for e^ = 1.2 e, e^ = e and e^ — 0.8 e. 
For each, we first considered the planar wall with poten- 
tial Vt^^z) = 2a^/z^ and determined the density profile 
PT^iz) and surface tensions ^wg, Iwi and 7 of the wall-gas, 
wall-liquid and liquid-gas interfaces, respectively. From 
Young's equation cos^ = (7^,^ — ^wi)/l we determined 
0{T) for each of these systems (see Fig. IT]). Each sys- 
tem exhibits a wetting transition, with T^, determined 
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FIG. 1: Variation of the contact angle with T for different wall 
strengths. The intersection with the dashed line at — 45° 
is the thermodynamic prediction for Tf. 



from the crossing of 7^,^ and jwi + 7. These occur at 
Tuj = 0.83 Tc, Tyj = 0.93 Tc and T^ = 0.99 T^ as e^ is 
reduced. The wetting transitions are all first-order; that 
is the thickness of the liquid layer i,^ jumps from a micro- 
scopic to macroscopic value at T^, . This is to be expected 
since the wall-fiuid potential is long-ranged but the trun- 
cated LJ fluid-fiuid interaction is effectively short-ranged 
2] . This prohibits critical wetting, which is important for 
our study. Also, as expected, the strength of the first- 
order transition decreases as T^, approaches Tc. This is 
apparent when one determines the interfacial binding in- 
terfacial binding potential W{£) corresponding to the ex- 
cess grand potential of a wetting film constrained to be of 
thickness £. The global minimum of this determines the 
equilibrium film thickness ^^. This is shown in Fig. [2] for 
the case e^ = 0.8 e close to the wetting temperature and 
shows an activation barrier, characteristic of first-order 
wetting, at £3 ~ 10 cr. For comparison the barrier for 
the binding potential for e^u = e is an order of magnitude 
larger and located at £b « 4cr. 

According to the thermodynamic prediction ([l]), the 
location of the filling transitions can be determined from 
the intersection of the contact angle curves with tp — 7r/4 
and gives T/ = 0.76 T^, Tj = 0.90 T^ and Tf = 0.97 T^ as 
Sw decreases in strength. To check this we set p — p~^^ 
and minimize 0[/5] to a global or local minimum fi, start- 
ing from different high density and low density configu- 
rations. For first-order filling these will converge to dif- 
ferent equilibrium profiles, corresponding to microscopic 
and macroscopic adsorptions, which coexist at Tf. This 
is what is found for the two strongest walls as illus- 
trated in Fig[3] where we plot the excess grand-potential 
fl°^ = i7 + pV per unit volume as a function of T. The 
values for Tf obtained are in near exact agreement with 
the thermodynamic predictions and differ from them only 
due to the limitations of numerical discretization and 
finite-size. In Fig. [4] we show the coexisting density pro- 
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FIG. 2; Binding potential function W{£) for e„ — 0.8 e close 
to a first-order wetting transition at T^ = 1.4 showing an 
activation barrier at Ib ~ 10 a. In the inset we show the 
binding potential at a lower temperature close to TJ ~ 1.38 
for which the activation barrier is still present. In both cases 
the results correspond to a bulk coexistence. 



files p{x,z), corresponding to macroscopic (left) and mi- 
croscopic (right) states, for e^ — £■ From these we can 
determine the thickness (.^ of the meniscus above the 
wedge apex defined as the distance from the origin to a 
point on a diagonal where p(x,x) = {pi + Pg)l^- The 
macroscopic meniscus is nearly flat (as it should be since 
we are at bulk coexistence) and meets each wall at the 
correct contact angle Q ss 7r/4. Of course the size of 
this macroscopic state is limited by our numerical grid 
and scales with the system size L. For the microscopic 
configuration the meniscus thickness i^ is larger than 
the wetting layer thickness £^ but of the same order as 
the distance of the activation barrier IIb ~ 4cr for the 
corresponding binding potential for the wetting transi- 
tion. This is precisely the expectation for first-order fill- 
ing from effective Hamiltonian theory [18]. Both micro- 
scopic and macroscopic profiles show layering behaviour 
close to the apex. 

For the weakest wall strength e^ = 0.8 e however both 
initial coverages converge to a unique phase indicating 
that the transition, which is of course rounded by the 
finite-size of our system, is continuous. A plot of the ad- 
sorption r — J J dxdz{p(x,z) — pg) versus T is shown 
in Fig. [5] and shows a dramatic but continuous increase 
in the adsorption near the anticipated Ti « 1.38. A 
cross-section of the density profile along the diagonal, 
p{x,x), for T ^ Tf is also shown (Fig. [6]). This indi- 
cates that the order of the filling transition is changed 
near the vicinity of the T^. Support for this comes from 
the two sources. Firstly near Tf the meniscus height 
£w « 22(7 is considerably larger than the location of the 
activation barrier is ~ 10 tr associated with the wet- 
ting binding potential (see inset Fig. [2]). This is not at 
all expected for first-order filling [lS\. Secondly we can 
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FIG. 3: Location of a first-order filling transition for Em ~ £■ 
Here V is the available volume which is the length of the 
wedge multiplied by {L — a)^ . 
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FIG. 4: Coexisting macroscopic (left) and microscopic (right) 
density profiles at a first-order filling transition for £„ — e. 



compare quantitatively with predictions for critical fill- 
ing. If the transition is critical then in an infinite wedge 
we expect that i^ ^ {Tf — T)^^™ with F ex £f^ owing to 
the triangular shape of the meniscus. Effective Hamilto- 
nian theory predicts that the critical singularities depend 
on the power-law describing the dominant wall-fluid or 
fluid-fluid interaction which we may write more gener- 
ally as V{z) « l/zP+^. The critical behaviour falls into 
two regimes with /?„, = 1/p for p < 4 and /S^, = 1/4 for 
p > 4 [TB]. Thus we anticipate /?„, = 1/2 in our model 
since p = 2. The inset in Fig. [5] shows a log plot of the 
adsorption for T < Tf, in which we use an unfitted esti- 
mate of the filling temperature TJ = 1.38 obtained from 
(fTl). This gives (3^ = 0.46 ± 0.05 in good agreement with 
the predicted value. 

The presence of a critical (or at least effectively crit- 
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FIG. 5: Temperature dependence of the adsorption in the 
wedge with the weakest wall interaction e,„ = 0.8 e. The inset 
shows the log-log plot of the adsorption vs the scaling field 
Tf - T. The slope of the straight line is -0.92. 




laboratory studies of complete filling in linear wedges 
corresponding to the approach to co-existence when the 
walls are completely wet {9 = 0) (5S] . However the obser- 
vation of critical filling would be more interesting because 
fluctuation effects arc much stronger. For example, we 
expect that beyond mean-fleld level, capillary-wave fluc- 
tuations do not alter the divergence of £^ « (T^ — T)~^/^ 
but do lead to a universal interfacial roughness (width) 
^j^ ~ {Tf—T)^^/^ which is much bigger than for complete 
fllling and also critical wetting [18]. These fluctuations 
are not captured by the present DFT and consequently 
the density profiles will be broader than predicted here. 
However our DFT should be otherwise extremely accu- 
rate regarding the location of the transition, its order and 
the adsorption. The observed change in order from first- 
order to continuous filling has only been partially antici- 
pated by previous effective Hamiltonian theory. This had 
been predicted on the basis of a simple interfacial model, 
valid only for shallow wedges, but the proposed mech- 
anism required that both the wall-fiuid and fluid-fluid 
forces to be of the same range. Then it was noted that 
even for flrst order wetting, the filling transition would be 
continuous if it occurs at a temperature below which the 
activation barrier forms in the binding potential W{i) 
[18j . However in the present DFT study a small activa- 
tion barrier is still present at Tj (see inset Fig. [2]), indi- 
cating that the prediction of the simple, shallow wedge, 
effective Hamiltonian theory is not completely correct. 
Nevertheless we believe that the substantial reduction in 
the size of the barrier as T approaches T^ plays a promi- 
nent role in the change of order of the filling transition. 
Finally it would be interesting to know if the change in 
order occurs via a tricritical or critical end point and also 
what happens for more acute wedges with stronger wall 
potentials. 



FIG. 6: Density profile along the diagonal from (0, 0) to (L, L) 
at the filling temperature Tf for e^ ~ 0.8 e. 



ical) filling transition when Tf is close to Tc, when the 
walls still exhibit first-order wetting, and in the presence 
of realistic long-ranged interactions is the main new result 
of our paper, and is we believe encouraging for experi- 
mental studies. To date there have only been detailed 



In this paper we have presented our results of numer- 
ical studies of first-order and critical filling transitions 
in a rectangular wedge using a non-local density func- 
tional theory. This is the first time that filling transitions 
have been studied using modern microscopic DFT in the 
presence of long-ranged wall-fluid interactions, and show 
that close to the bulk critical temperature the wedge fill- 
ing transition may be continuous even though the walls 
themselves exhibit first-order wetting. 
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